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We present the bulk thermodynamic properties and phase diagram of strongly interacting matter 
in an extension of the 3-flavor NJL and PNJL models of QCD. Using a three momentum cut-off 
scheme, we have extended the multiquark interaction terms up to eight order so that the stability 
• of the vacuum is ensured in these models. We explore the effects of various combinations of the 

Oh' two eight-quark couplings gi and g2 and present a comparative study between the NJL and PNJL 

models as well as Lattice QCD data. The main effect of the eight-quark interaction term is to shift 
the critical end point in the T — /i phase diagram to a lower value of /i and higher value of T, thus 



(— I ' bringing them closer to Lattice QCD results. 

PACS numbers: 12.38.Aw, 12.38.Mh, 12.39.-x 
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■ I. INTRODUCTION 

rn 

CO ' The strong interaction as described by Quantum Chromodynamics (QCD), is a remarkable 



branch of physics which shows a rich phase structure at finite temperature and density. At 
low temperatures and densities, the dominant degrees of freedom in our nature are color-singlet 
bound states of hadrons. However, due to the asymptotic freedom of QCD, it is expected that at 
very high temperatures and densities these hadrons break up to liberate quarks and gluons and 
k>( , form the quark-gluon plasma (QGP). The experimental exploration of such a phase transition 

j_j ■ from the confined hadronic phase to the deconfined QGP phase is being pursued actively in 

' the Rclativistic Heavy Ion Collider (RHIC) and more data are expected from Large Hadron 

Collider (LHC) running at CERN and in the future experiment at FAIR. 

The transition region is however quite far from the asymptotic regime of QCD. This non- 
perturbative nature makes the study of these hot and dense matter quite non-trivial. The 
most reliable way to analyze the physics in this range of interest is to perform the numerical 
computation of a lattice version of the color SU(3) gauge theory (Lattice QCD/LQCD). The 
scheme is robust but numerically costly. Hence the most popular way to study the physics of 
the phase transition of strongly interacting matter is to look into one of the various possible 
effective models inspired from QCD. One of these is the Nambu-Jona-Lasinio (NJL) model [ll-Q 
which incorporates the global symmetries of QCD quite nicely. A four quark interaction term in 
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the NJL Lagrangian is able to generate the physics of spontaneous breaking of chiral symmetry 
- a property of QCD which is manifested as the nondegenerate chiral partners of the low-mass 
hadrons. But the major drawback of the NJL model is to have a reasonable description of the 
physics of color confinement. In this respect, the Polyakov loop extended NJL (PNJL) model 
ja-Hl tries to incorporate the fact that the chiral transition in QCD is of deconfining nature by 
introducing a background temporal gluon field. 

In our previous work [9|] we elaborately studied the 2+1 flavor PNJL model [Si], incorporating 
four-quark and six-quark interaction terms in the Lagrangian with three-momentum cut-off 
regularization. The six-quark interaction term (or the 't Hooft determinant term) preserves 
the SU{3)j^ X SU{3)j^ symmetry and breaks the C/(l)^ symmetry required due to the chiral 
anomaly. This term is also responsible for the flavor mixing in the tjq and rjs mesons to give rj 
and 77' mesons. 

However, due to the six-quark term the vacuum of the model becomes unbounded from below. 
This problem has become more obvious in the functional analysis of the Lagrangian. In the 
work of Reinhardt and Alkofer [lol | the stationary phase approximation method was considered 
to bosonize the model and to calculate the contribution of the classical path at the lowest order. 
This lowest order results sums all tree diagrams in powers of the coupling constant of the six 
quark interaction. But the functional integral bosonization of the model shows several classical 
trajectories in that interval of the functional integration. Since at leading order there can be 
only one classical trajectory in the mean-field potential, the semiclassical potential derived from 
the functional integration differs from the mean-field one. So if we consider several classical 
trajectories, the effective potential of the theory gets unbounded from below. 

In fact, the origin of multi-quark interactions is not very well-known. But we can find 
evidences of multi-quark interactions in the semi-classical theories based on QCD instanton 
vacuum [flj . There we can find certain correlations between two or more quarks by averaging 
over their positions and orientations in color space. Also in the instanton-gas model an infinite 
number of multi-quark interactions can be found beyond the zero mode approximation (l^ . 
There are also some lattice measurements for the QCD vacuum which shows a hierarchy between 
the multi-quark interaction |13ij . In this case the lowest four quark interaction term forms a 
stable vacuum by breaking the chiral symmetry spontaneously. But the next term in the 
hierarchy, the six quark interaction term, which is needed to mimic the J7a(1) anomaly, destroys 
the ground state. So one cannot truncate the tower of multi-quark interactions at this level. The 
next candidate is the eight-quark interaction term which ensures the stability of the vacuum. 

There are some recent findings [l3 - [T6| which show that the addition of the eight-quark 
interaction term in the Lagrangian may solve the problem of unstable vacuum. The (2-1-1)- 
flavor NJL model with an eight-quark interaction was studied in the four momentum cut-off 
and Pauli-Villars regularization scheme. The standard bosonization procedure of six and eight- 
quark interactions was followed and the multi-quark vertices were replaced by purely mesonic 
ones by the stationary-phase method. Indeed the controversy between the mean-field approach 
and the functional integral approach has been removed by including the eight-quark interaction 
term in the Lagrangian, since it restricts the number of classical trajectories to one, giving a 
stable ground state of the system. Recently, Kashiwa et. al. [l3, [13 have studied the chiral 
phase transition in the 2-fiavor NJL and PNJL models with eight-quark interactions. The 
2-flavor model does not incorporate the six-quark interaction. 

In this paper we first of all extend the (2-|-l)-flavor NJL and PNJL models to include the eight- 
quark interactions at non-zero temperatures and densities. For this we use the 3-momentum 
cut-off regularization scheme. Thereafter we thoroughly investigate the effect of the eight-quark 
term in the thermodynamic properties of strongly interacting matter. 

We organize our paper as follows: in the next section we briefiy describe the 2-1-1 fiavor PNJL 



3 



model with eight-quark interaction term in the Lagrangian. For the Polyakov loop potential 
we incorporate Vandermonde term In sec. /// we describe the thermodynamic potential of 
the PNJL model. We elaborately discuss the parameter sets chosen, and the stability criteria. 
The next two sections discuss the results of the paper, about the chiral phase transition and the 
bulk thermodynamic properties like pressure, trace anomaly at zero density and quark number 
density at finite temperature and density. We give a detailed analysis of the phase diagram and 
the critical end points at the next section. A summary of our results are available in the last 
section. 



II. THREE FLAVOR PNJL MODEL WITH EIGHT-QUARK INTERACTION 

In this section wc briefly describe the formalism of the PNJL model with eight-quark inter- 
action term. Some details on three flavor PNJL model with four and six quark interaction can 
be found in the literature [1, @, d, and the NJL model with eight-quark interaction term in 
There is also some work on two flavor PNJL model with eight-quark interaction term 



In the PNJL model the gluon dynamics is reduced to the chiral point couplings between quarks 
and a simple temporal background gauge field which represents Polyakov Loop dynamics. The 
Polyakov line is represented as (see e.g. [l^ and references therein), 

= Pexp[i / dTAi{x,T)] (1) 



^0 

where = iAo is the temporal component of Eucledian gauge field (A, A4), /3 = y, and V 
denotes path ordering. L{x) transforms as a field with charge one under global Z(3) symmetry. 
The gluon dynamics can be described as an effective theory of the Polyakov loops 120|. The 
quark thermodynamics can be described by the NJL model Let us consider the 5'{/(3)/ 
version of PNJL model with eight-fermion interaction described by the Lagrangian, 

f=u,d,s f f 

0=0,. ..,8 

-W'($[A],I>[A],T) 

= £0 + ^SB + + Cs + CjCMT + + t^\q - ^' (2) 

where the matrices Pl,r — (1 i 75)/2 are chiral projectors. In the above Lagrangian is the 
Dirac term with gauge field interactions; = — iA^S ^4^. Csb is the mass term which breaks 
the symmetry explicitly. The mass of a particular flavor is denoted by my, where f — u,d, s. Cs 
is the four-fermi interaction term with coupling gs- The next term, Ckmt, which is a six-fermi 
interaction term, is invariant under SU{i)j^ x SU{i)p but breaks U{\)a symmetry. This term 
mimics the QCD chiral anomaly. The terms and are the eight-quark interaction terms 
which describe the spin zero interactions where gi and 52 are the corresponding couplings. Here 
we have considered the interaction terms to be effectively local since the meson physics in the 
large limit is described by the local Lagrangian of this type. Since the coupling constants 
are dimensionful, the model is not renormalizable. So we have used three-momentum cut-off 
regulator A to make quark loops finite. 
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In earlier versions of the PNJL model we found that the Polyakov loop, which is the nor- 
malized trace of the Wilson line L, has become greater than 1 above 2Tc [1, HH, [l^]- To solve 
this problem one has to take a proper Jacobian of transformation from the matrix valued field 
L to $ which will then constrain the value of $ within 1. One way to resolve this problem is 
to introduce Vandermonde term in the Polyakov loop potential [3] ■ Thus the potential W with 
the Vandermonde term can be expressed as 

^'('^[-^]^f[-^]'^) . - . ln[ J($, $)] (3) 

where U{(j)) is the Landau-Ginzburg type potential given by [Q], 

with 

$ = {Tr,L)lN,, $ = {Tr,L^)/N, 
b^{T) = ao + ai(^) + a^i^f + a,{^f , (5) 

63 and 64 being constants. The second term in eqn. |3] is known as Vandermonde term, where 
J($, is the Jacobian of transformation from Wilson line L to ($, i>) written as 

J[$, !>] = (27/247r2)(l - 6$l> + 4($^ + 1-^) - 3($$)^ 

J($, is also known as Vandermonde determinant and is not explicitly space time dependent. 
The value of the dimensionless parameter k will be determined phenomenologically. The co- 
efficient K in the VDM term can in general have some temperature and/or chemical potential 
dependence. But here we take a constant value in such a way that we can get the pressure and 
the transition temperature as close as possible to the lattice QCD results for the PNJL model 
with and without the eight-quark interaction. This coefficient has been tabulated later. 

In order to study the chiral transition of the system we consider the Mean Field Approx- 
imation (MFA) of the eqn. ([2]) to get the field equations for $, cr, where it is assumed 
that the system is described as an assembly of non-interacting particles moving in the mean 
field. The theory is analogous to the BCS theory of superconductor, where the pairing of two 
electrons leads to the condensation causing the gap in the energy spectra. Similarly in this 
model, due to the dynamical breaking of S'C/(3)^ x SU{?>)r symmetry to SU{?>)v a composite 
operator picks up a nonzero vacuum expectation value leading to (V'V') condensation. Due to 
the dynamical breaking of chiral symmetry, nine goldstone bosons appear for Nf = 3 model. 
The quark condensate is given as, 

= -iNcCty^^+{trSf{x - y)) (6) 

where trace is over color and spin states. The self-consistent gap equation for the constituent 
masses are, 



Mf =mf - 2gsaf + — cr/+if7/+2 - 2giaf{al + + cr^) - 452^^ (7) 



here cr/ = (tl'f'fpf) denotes chiral condensate of the quark with flavor /. / denotes the flavor 
u,d,s respectively. Here if we consider af — cr„, then (Jf+i = (Jd and iyf+2 = Cs- Similarly 
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if CT/ = Gd then (Jf+i = Os and (Jf+2 — <^d, if cr/ = cts then cr/+i = (Ju and cr/+2 = cTd- The 
expression for af aX T — Q and ^ — Q can be written as [8] 



3M, 



A being the three-momentum cut off. 



-.dp, 



(8) 



III. THERMODYNAMIC POTENTIAL AND PARAMETER FITTING 



The primary aim of our work is to study the thermodynamic properties of the strongly inter- 
acting matter using the PNJL model with eight-quark interactions at zero and finite chemical 
potentials. To do so we now need to fix the parameters in both the NJL and Polyakov loop 
potentials. The thermodynamic potential for the multi-fermion interaction in MFA of the PNJL 
model can be written as 



r!($,$,M,r,^) =W'[$,l>,r] + 2.gs ~ '^<^u<Tdas+i^{(T)f + ig2(T^f 



(27r)3 



2 2 

f—u.d,s 



T 



(9) 



where a;„ = 7rT(2n+ 1) are Matsubara frequencies for fermions. The inverse quark propagator 
is given in momentum space by 



using the identity Tr In {X) = Indct (AT), we get 



(10) 



f—u^d,s 
, Jo (2TT) 



^ (2^) 



-f Jo i2ny 
<Pp 



2^T 



7^ Jo i2nr 



■In 



In 



1 + 3($ + $e r )e T -i-e r 



l + 3($-|-$e T 'je T -(- e t 



(11) 



where Epf — y^p2 + M| is the single quasiparticle energy, — (ct^ + <^d + '^s) ^^d aj — 

(cr^ + a'^ + aj). In the above integrals, the vacuum integral has a cutoff A whereas the medium 
dependent integrals have been extended to infinity. The main idea to incorporate the eight- 
quark interaction term in the Lagrangian is to stabilize the vacuum. The detailed discussion of 
the stability criteria of the vacuum can be found in the literature [IBl • 

The parameters of the NJL part of the Lagrangian are the current quark masses to«, rrid and 
TTis , the coupling constants gs, go, gi and 172 and the three- momentum cutoff A characterizing 
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Sets 


Physical 
Parameter 
values 


six-quark 


eight-quark 


a 


b 


c 


d 


e 


f 


g 


h 


















= 138 MeV 


U 






U 


U 






f. 


= 93 MeV 


ruk 


rrik 


ruk 


rrik 




mu 


rrik 


rrik 


rUk = 494 MeV 










fk 


h 


fk 


fk 


fk = 117 MeV 


rrin 


rrir, 


rrin 




nin 




rrir, 




= 480 MeV 










m^, 






m^, 


m^, = 957 MeV 




rria 














m„ = 680 MeV 








ruu 


ruu 






ruu 


rriu = 5.5 MeV 














(Ju 




cry^ = -248 MeV 



TABLE I: (Color online) The PNJL model parameters are fitted to reproduce the different physical 
parameters as indicated above. 



Set 


mu 
MeV 


MeV 


A 
MeV 




goA-" 


gi X 10"^^ 
MeV"* 


52 X 10"^^ 
MeV"** 




rpPNJL 

MeV 


rpNJL 

MeV 


a 


5.5 


134.758 


631.357 


3.664 


74.636 


0.0 


0.0 


.13 


181.0 


170.55 


b 


5.406 


133.227 


641.357 


3.717 


61.309 


0.0 


0.0 


.18 


182.2 


170.25 


c 


5.418 


133.562 


640.206 


3.637 


70.849 


0.0 


0.0 


.14 


180.5 


169.25 


d 


5.5 


133.532 


631.337 


4.229 


14.461 


0.0 


0.0 


.11 


183.5 


176.65 


e 


5.5 


183.468 


637.720 


2.914 


75.968 


2.193 


-5.890 


.06 


168.5 


141.35 


/ 


12.509 


181.863 


628.933 


2.986 


75.444 


2.007 


-4.538 


.07 


171.8 


150.05 


9 


8.742 


179.498 


640.206 


2.928 


75.382 


1.929 


-4.840 


.05 


169.0 


144.25 


h 


5.5 


187.786 


628.933 


2.956 


75.983 


2.425 


-6.445 


.08 


170.5 


144.15 



TABLE II: (Color online) Parameters of the 5(7(3) NJL part and different value of Tc at ^ = and 
the value of k for all input parameter sets. 



the scale of the chiral symmetry breaking. Here we consider the isospin symmetric limit m„ = 
TTid. These parameters are determined by reproducing few physical quantities like pion mass 
jTiTT = 138 MeV, kaon mass uik = 494 MeV, eta mass to,, = 480 MeV, eta prime mass to^ = 
957 MeV, pion decay constant /^r = 93 MeV, kaon decay constant Jk = 117 MeV, u quark 

1 /3 

condensate cr„ = —248 MeV and/or the mass of sigma meson TOo- = 680 MeV. 

We consider two cases of up to six-quark (gi = = 32) and up to eight-quark (gi ^ 7^ 172) 
interactions each with four sets (a, 6,c, d) and {e,f,g,h) respectively. In table H] we have tab- 
ulated different physical observables and their values for the different sets. In set (a, d) and 
(e, h) rriu is kept fixed at 5.5 MeV at 1 GeV scale since the values determined from the eval- 
uation of the current matrix elements at low energies are also centered around 5.5 MeV. In 
sets b and / the mass of sigma meson is considered as fixed to obtain the model parameters 
whereas in sets c and g the u quark condensate is considered. In sets (b, c) we fit the parameters 
(mujiris, A, gs, go) by fixing (to^, m/f, m,,, /tt and nia- or tT„) and in sets (/,(?) the parameters 
(m„,TOs, A,gs,g_D,5i,52) are obtained by fixing {m^^mK^mj^.m' f^^^K and m^r or ct„). How- 
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ever in set d and h instead of we have used mg- — 680 MeV for fitting tlie same parameters 
as Too- is the most sensitive to the eight-quark couphngs. 

The sensitivity of sigma mass to the eight-quark couphng constant have been also discussed 
in the paper [15j, however instead of producing the eight-quark couphngs, they put the values 
of these coupling constants by hand. The reason behind producing four sets of parameters is 
to show the effect of different physical channels on the different parameters. We would like to 
point out here that the mass of the sigma meson is taken as 680 MeV to fit the parameters 
in sets b and /, due to which we obtain m„ — 12.5 MeV in set /. The value of mu in set / 
is higher than the value of m„ = 5 — 9 MeV that is usually quoted in the literature. A lower 
mass value of sigma meson ^ 660 MeV may give a lower value of current u quark mass but 
then the constituent u quark mass will be much higher than the value obtained for other input 
parameter sets. In our analysis we will consider the plots for the input parameter sets in the 
following combinations: 

• Set 1 : Sets (a, e) (six-quark, eight-quark) . Here — 5.5 MeV is held fixed. 

• Set 2 : Sets {b, f) (six-quark, eight-quark) . Here rrio- = 680 MeV is used to fit parameters. 

• Set 3 : Sets {c,g) (six-quark, eight-quark) . Here (T„ is used instead of mg-. 

• Set 4 : Sets {d,h) (six-quark, eight-quark) . Here m„ — 5.5 MeV is held fixed and rria — 
680 MeV is used to fit parameters. 

The resulting model parameter values of the NJL part at zero temperature, obtained in our fit 
is given in the columns 2-8 of table HIl 

For the Polyakov looppart we have to consider the finite temperature behavior of the PNJL 
model. As discussed in [7] we are able to tune the dimensionless coupling k in the Vandermonde 
term and obtain reasonable behavior of the mean fields. We thus choose the following set of 
parameters, 

ao = 6.75, ai = -1.95, 02 = 2.625a3 = -7.44, 63 = 0.75, 64 = 7.5, To = 190 MeV 

The remaining parameter is the dimensionless coupling k in the Vandermonde term which has 
been tabulated in column 9 of table [Hi This has been obtained by choosing suitable value of k 
so that the pressure in PNJL model follows the Lattice QCD data as closely as possible. 

Our job is now to estimate the transition temperatures at different quark chemical potentials. 
In order to obtain this we search for the minimum of the thermodynamic potential which gives 
the temperature and density dependence of the fields. The point of inflection or the gap, as 
the case may be, of these fields give the transition temperature and chemical potential. This 
tusk has been carried out for both NJL and PNJL models. The transition temperatures at zero 
chemical potential have been tabulated in the last two columns of table |IT1 The most important 
observation in this regard is that the introduction of the eight quark interaction lowers the 
transition temperature to the range 150 — 190 MeV. 

The pressure of the strongly interacting matter obtained from our expression of the thermo- 
dynamic potential is, 

P{T,fi)^-n{T,pL), (12) 

where T is the temperature and /x is the quark chemical potential. In fig. [T] we show the 
variation of the scaled pressure P/Psb with T/Tq in both the NJL and PNJL models for up to 
six-quark as well as up to eight-quark interactions. The plots show that the pressure is slightly 
smaller near the transition region due to eight-quark interaction term. In general, the shift of 
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pressure in the whole range of temperature raiiges from 10% — 20%. We have also compared 
our results with those of the lattice QCD data [23| with temporal extent Nt = 4 and Nt — 6. 
In fact the pressure in the PNJL model was obtained to give the best possible fit to the lattice 
data with Nt = 6. We find that for all input parameter sets it is possible to have an impressive 
similarity with the results of the lattice data. On the other hand the results for the NJL model 
differ considerably from the lattice data. 

Another basic thermodynamic quantity is the energy-momentum tensor. The thermal con- 
tribution of the trace of this energy-momentum tensor can be defined by the difference of e and 
3p, where e is the energy density and p is the pressure. So in thermodynamics the traceless 
energy-momentum tensor means = (e — 3p) ~ 0. The trace of the energy momentum 
tensor is vanishing at the classical level when the theory has no mass scale. We know that 
QCD in the chiral limit is scale invariant, which means that in massless QCD, is zero 
unless quantum corrections are taken into account. Thus in a conformally symmetric theory, 
for example a theory of free massless gluons, = 0. Therefore this quantity measures the 
breaking of conformal symmetry in the interaction theory. It is straightforward to evaluate the 
(e — 3p) or the interaction measure from the thermodynamic potential and pressure as follows: 

(. - 3p)/r^ = T^-^^ (13) 



In fig. [5] we show the variation of the (e — 3p)/T^ with temperature for both the models 
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FIG. 2: (Color online) Variation of (e - 3p)/r* with T/Tc, at /i = for both NJL and PNJL model. 
For details of the sets see text in section [Till 



with and without eight-quark interaction. At temperatures just above Tq the plots show a 
peak due to the largest deviations of (e — 3p) from the conformal limit, e — 3p. These peaks 
establish a prominent structure of trace anomaly which is consistent with lattice calculations. 
The introduction of finite value of gi and 172 slightly lowers the peak position in PNJL model 
which ensures reduction of trace anomaly due to the eight-quark interaction term. At high 
temperatures the trace anomaly dropped rapidly similar to that obtained from the lattice 
data. From the plots we can see that in this temperature region much larger than Tc, the 
trace anomaly drop less rapidly for plots with finite value of gi and 172 than the plots with 
51 = 32 = 0. In all diagrams the peak position for Nr = 4 result is much higher than the 
Nr = 6 and the peak positions of the (e — 3p) /T"* for the PNJL model are slightly lower than 
the peak position of the lattice result for Nr = 6. It remains to be seen if this is close to the 
continuum limit of lattice data. At low temperature region the plots for all the parameter sets 
coincide very well with the lattice data, but at high temperature the parameter sets with the 
six-quark interaction are in better agreement with the lattice data. However the deviation of 
the plots with the eight-quark interaction term from the lattice results are also very small and 
a continuum extrapolation can be expected to yield even a better agreement. In pure gauge 
theory, it has been noted that (e — 3p) is quite significant above the deconfinement temperature 
[23] and ~ for temperatures up to a few times the deconfinement temperature ^25i] . A similar 
behavior can be found in the PNJL model. 
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FIG. 3: (Color online) Variation of Uq/T'^ with temperature, at /i = 100 MeV for both NJL and PNJL 
models. For details of the sets see text in section UlTl 



The scaled quark number density is defined as: 

We have plotted the quark number density as a function of T/Tc a.t fi ~ 100 MeV for both 
the models in fig. [3] and studied the effect of eight-quark interaction. It can be seen that at 
fixed values of temperature and chemical potential the Uq in PNJL model is much lower than 
that in the NJL model below Tc. This is evidently an effect of confinement in the PNJL model. 
At a fixed quark chemical potential the quark number density for the PNJL model is almost 
vanishing below the chiral transition temperature, but rises very quickly in the vicinity of the 
transition temperature. However in case of the NJL model the rise in the quark number density 
starts much below the transition temperature. Interestingly, the addition of the eight-quark 
interaction term raises the quark number significantly above Tc for all cases at hand. 



IV. PHASE DIAGRAMS AND CRITICAL END POINT 



The position of the critical endpoint (CEP) in the phase diagram is one of the most inter- 
esting issue of the hot and dense strongly interacting matter. The phase diagrams are usually 
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obtained by identifying the critical temperature with the temperature at which the light quark 
chiral condensate has a jump (first order) or from the maximal point of the derivative of light 
quark condensate with respect to the temperature for different chemical potentials. CEP is 
the point which separate the cross-over transition from the first order phase transition. There 
have been a number of studies in lattice QCD to find out the CEP in the T — fj, diagram 
[26*1 and also in QCD inspired models. In three flavor NJL model with the parameters of 
Hatsuda-Kunihiro 3] the location of CEP is found to be (Tcfic) = (48MeV, 324 MeV), how- 
ever in our previous work with three flavor PNJL model we have shown the position of CEP 
at {Tc,^ic) = (92MeV, 314MeV) which establishes the fact that the location of the CEP goes 
up in temperature in the PNJL model The reason behind this shift to the higher value of 
temperature is due to the suppression of quark excitations at finite temperature and density by 
the Polyakov loop. However the lattice estimates of the CEP vary cosiderably between different 
groups and the value of /ic (< 150 MeV) tend to be much lower than the PNJL value. In this 
paper we try to develop a comparative study with our four input parameter sets for both NJL 
and PNJL model with and without eight-quark interaction. In fig. 2] and fig. [S]we have shown 
the phase diagrams of NJL and PNJL model with and without eight-quark interaction term for 
our four input parameter sets. In the region T less than Tc and /i greater than fj,c the chiral 
and deconfinement transitions are first order and occur almost at same T and ^. This can be 
realized from the generalized Clausius-Clapeyron relation for the system with multiple order 
parameters, which shows that in case of first order phase transition all discontinuities appear 
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FIG. 5: (Color online) Phase diagram in /i with T for PNJL model. For details of the sets see text in 
section IIIII 

at same T and /i. In table IIIll we have shown the values of CEP for different sets of parameters 
in both NJL and PNJL model. From the parameter sets we can clearly infer that the newly 
added eight-quark interaction term shifts the CEP towards lower fi and higher T value for 
both NJL and PNJL model. The recent work in lattice QCD predicts that the possible region 
where CEP should exist is fic/Tc < 2.5 From the table Hill we can observe that the PNJL 
model with eight-quark interaction term only satisfy this criteria. For the PNJL model without 
eight-quark interaction term this ratio becomes greater than 3 for all parameter sets. Moreover 
in case of NJL model with eight-quark interaction term, this ratio is not within the allowed 
region. This feature establishes the importance of including the eight-quark interaction term 
in the Lagrangian. 

V. CONCLUSION 

To conclude, we have studied the bulk thermodynamic properties and the phase diagram 
of the PNJL and NJL model with eight-quark interaction. This is the first time that 3-flavor 
PNJL model has been studied with eight-quark interaction and three-momentum cutoff. 

We have used different sets of physical observables to fix the input parameters for our model 
(including the eight quark couplings). In particular, we can observe the sensitivity of the sigma 
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Model Index 


inc MeV, Tc MeV) 


set 1 


set 2 


set 3 


set 4 


NJL (6 - quark) 
NJL (8 - quark) 
PNJL (6 - quark) 
PNJL (8 - quark) 


(323,48.5) 
(263,61.2) 
(313,92.85) 
(260, 118.5) 


(325,40.15) 
(305,46.55) 
(313,88.5) 
(237, 122.05) 


(325,39.8) 
(283,51.4) 
(317, 82.25) 
(255, 121.95) 


(328, 46.45) 
(263,65.95) 
(320,89.15) 
(213, 134.75) 



TABLE III: (Color online) The values of (AtcTc) for different sets of NJL and PNJL model 



mass to the eight-quark coupling constants. However the value of sigma mass mo- = 680 MeV 
we took as one of the variable for fitting the input parameters with eight-quark interaction term 
produce niu — 12.5 MeV which is higher than the range 5 — 9 MeV. 

The effect of eight-quark interaction on thermodynamic properties are the main interest of our 
work. The variation of scaled pressure P/Psb with T /Tc for four input parameters have been 
plotted and the results are compared with available lattice data. The inclusion of eight quark 
interaction decreases the pressure near the transition temperature and increases the pressure 
at higher temperature. We have also studied the variation of (e — 3p)/T'^ with temperature. 
The six-quark interaction is already in good agreement with N^- — 6 lattice data. However, 
comparison with N-r = 4 lattice data gives us the impression that a continuum extrapolation 
may give a better match to our results for eight-quark interactions. 

We discuss the effect of higher order interaction term on the finite density results. We observe 
that due to the eight-quark interaction quark number density increases above the transition 
temperature. 

The effect of the eight-quark interaction term on the critical end points in the phase diagram 
is the main issue of our study. From table IHII we can observe that the eight-quark interaction 
drives the CEP to a low chemical potential and a high temperature value. Furthermore, our 
study concludes that the inclusion of eight-quark interaction is essential to limit /ip /Tq below 
2.5 as suggested by Lattice calculations. 

All these studies lead us to the conclusion that the eight-quark interaction term in the La- 
grangian has an interesting phenomenological implication in effective models. 
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